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Time-series  measurements  of  daylight  (skylight  plus  direct  sunlight)  spectra  beneath  overcast  skies 
reveal  an  unexpectedly  wide  gamut  of  pastel  colors.  Analyses  of  these  spectra  indicate  that  at  visible 
wavelengths,  overcasts  are  far  from  spectrally  neutral  transmitters  of  the  daylight  incident  on  their  tops. 
Colorimetric  analyses  show  that  overcasts  make  daylight  bluer  and  that  the  amount  of  bluing  increases 
with  cloud  optical  depth.  Simulations  using  the  radiative-transfer  model  MODTRAN4  help  explain  the 
observed  bluing:  multiple  scattering  within  optically  thick  clouds  greatly  enhances  spectrally  selective 
absorption  by  water  droplets.  However,  other  factors  affecting  overcast  colors  seen  from  below  range  from 
minimal  (cloud-top  heights)  to  moot  (surface  colors).  ©  2005  Optical  Society  of  America 
OCIS  codes:  010.1290,  290.1090,  330.1710,  330.1730. 


1.  Introduction 

In  the  extended  family  of  atmospheric  optical  phe¬ 
nomena,  overcast  skies  are  surely  the  poor  relations. 
Lacking  the  distinct,  vivid  colors  of  such  prodigals  as 
parhelia  and  rainbows,  overcasts  have  prompted 
relatively  little  spectral  research  outside  the  ultravi¬ 
olet  and  infrared,  where  their  radiative-transfer 
properties  are  important.1-3  With  relatively  few 
exceptions,4-6  even  those  researchers  who  do  analyze 
overcasts  at  visible  wavelengths  are  interested 
chiefly  in  explaining  the  clouds’  spectral  properties 
independent  of  human  perception.3-7-9  And  although 
the  CIE  has  recommended  a  particular  model  for  the 
angular  distribution  of  overcast  luminances,10  it  is 
far  less  specific  about  the  spectral  distribution  of 
overcast  radiances  or  irradiances,  instead  indicating 
that  overcast  spectra  are  simple  variants  on  its  stan¬ 
dard  daylight  illuminant  D65. 1 1 

This  relative  scientific  neglect  of  overcast  color  is 
paired  with  popular  disdain:  in  everyday  parlance, 
overcasts  are  the  dark,  dreary,  and  dull-colored  al¬ 
ternatives  to  visually  pleasing  clear  or  partly  cloudy 
skies.  Implicit  in  such  language  is  the  plausible  (but 
incorrect)  assumption  that  overcasts  are  spectrally 
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neutral  transmitters  of  the  daylight  (skylight  plus 
direct  sunlight)  incident  on  them.12  Yet  our  research 
shows  that  overcast  colors  are  seldom  as  gray  as  we 
might  imagine,  and  explaining  why  overcasts  are 
spectrally  selective  yields  useful  insights  into  radia¬ 
tive  transfer  in  clouds. 

2.  Measured  Chromaticities  of  Daytime  Overcasts 

Because  “overcast”  has  slightly  different  meanings  in 
different  contexts,  we  start  by  defining  the  term.  In 
our  work,  an  overcast  must  meet  two  criteria:  ( 1)  no 
clear  sky  can  be  visible  anywhere  and  (2)  cloud  cover 
must  be  sufficiently  optically  thick  that  any  cast 
shadows  are  indistinct.  Unlike  some  definitions  of 
overcast,13  ours  includes  surface  fog  that  obscures  (or 
perhaps  is)  the  sky.  We  restrict  ourselves  here  to 
daytime  overcasts  when  unrefracted  Sun  elevation 
h0  >  5°  because  overcast  color  behaves  quite  differ¬ 
ently  at  lower  h0.  Although  most  of  our  overcasts  were 
formed  by  thick  stratus  (St)  or  stratocumulus  (Sc) 
clouds,  we  also  measured  altostratus  overcasts 
through  which  the  Sun’s  disk  was  faintly  visible.  To 
protect  our  instruments,  we  measured  mostly  non¬ 
precipitating  overcasts.  With  care,  however,  we  also 
were  able  to  acquire  some  data  in  drizzle,  light  rain, 
and  snow.  What  kinds  of  chromaticity  gamuts  do 
such  overcasts  produce? 

Figure  1  shows  part  of  the  CIE  1976  uniform- 
chromaticity-scale  (UCS)  diagram,14  within  which  we 
plot  temporal  chromaticity  trends  produced  by  two 
typical  overcasts.  We  generate  such  curves  by  mea¬ 
suring  overcast  spectra  every  30  s  and  then  connect¬ 
ing  in  sequence  the  individual  u',v'  chromaticity 
coordinates  calculated  from  each  spectrum.  Included 
as  colorimetric  landmarks  in  Fig.  1  are  part  of  the 
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Fig.  1.  Portion  of  the  CIE  1976  UCS  diagram,  showing  temporal 
trends  in  stratus  (St)  and  stratocumulus  (Sc)  overcast  chromatici- 
ties  measured  at  the  U.S.  Naval  Academy  in  Annapolis,  Maryland 
(USNA)  on  4-21-2003  and  at  the  University  of  Granada  in 
Granada,  Spain  on  11-8-03.  Chromaticities  are  calculated  from 
horizontal  spectral  irradiances  Ex.  For  comparison,  we  also  show 
chromaticities  from  two  cloudless  days  at  USNA  (9-5-02,  9-10-03). 
On  each  day,  unrefracted  Sun  elevation  h0  ranges  from 
30.0°-36.3°. 


Planckian  locus  and  its  corresponding  color- 
temperature  limits.  Figure  l’s  observing  sites  are  the 
University  of  Granada  at  Granada,  Spain,  and  the 
U.S.  Naval  Academy  (USNA)  at  Annapolis,  Mary¬ 
land.  We  selected  the  Fig.  1  curves  from  our  time 
series  of  >9100  radiance  and  irradiance  spectra  mea- 


Table  1.  Geographic  Details  of  Our  Measurement  Sites 


Site  Name 

Latitude 

Longitude 

Elevation  (m) 

Granada,  Spain 

37°  11'  N 

3°  35'  W 

680 

Marion  Center,  Pa. 

40°  49'  N 

79°  5' W 

451 

Owings,  Md. 

38°  41'  N 

76°  35' W 

15 

USNA,  Md. 

38° 59' N 

76° 29' W 

18 

sured  during  40  overcasts  at  USNA  and  Granada,  as 
well  as  at  two  rural  sites  near  Owings,  Maryland,  and 
Marion  Center,  Pennsylvania.  Table  1  gives  details  of 
the  sites’  locations;  we  made  all  observations  at 
ground  or  rooftop  level. 

Our  instruments  are  two  Photo  Research  PR-650 
spectroradiometers,  with  which  we  measure  either 
zenith  spectral  radiances  Lh  or  horizontal  spectral 
irradiances  Ex  at  visible  wavelengths  X.  Both  radiom¬ 
eters  have  spectral  ranges  of  380-780  nm  and  a  step 
size  of  4  nm.15  Zenith  radiances  are  measured  across 
a  1°  field  of  view  (FOV)  by  using  a  telescopic  lens  on 
the  PR-650,  and  hemispheric  irradiances  are  mea¬ 
sured  using  a  nominally  cosine-corrected  diffuser 
(FOV  =  2'n  sr).  We  aimed  the  radiometer  and  then 
locked  it  on  a  tripod  before  each  measurement  ses¬ 
sion.  Figure  l’s  chromaticities,  are  calculated  from 
horizontal  Ex  and  show  only  parts  of  such  sessions: 
each  curve  spans  the  period  when  30°  <  h0  <  36.3°. 
Figure  1  also  shows  the  range  of  clear  daylight  chro¬ 
maticities  (abbreviated  CLR  and  plotted  using  x’s) 
measured  on  two  clear  days  for  the  same  h0  interval. 
The  conventional  meteorological  abbreviations 
“OVC”  and  “CLR”  in  Fig.  1  indicate  overcast  and 
clear,  respectively. 

Several  colorimetric  features  are  noteworthy  in 
Fig.  1.  First,  its  chromaticity  gamut  for  clear  daylight 
is  much  less  than  for  its  two  stratus  and  stratocumu¬ 
lus  overcasts.  Using  the  normalized  colorimetric 
gamut  g,m  Fig.  l’s  two  overcast  gamuts  are  7.0  and 
13.8  times  larger  than  the  combined  clear-sky  gamut 
(see  Table  2).  Thus  an  overcast’s  colors  are  only  partly 
determined  by  the  daylight  color  that  illuminates  its 
cloud  tops.  Second,  Fig.  1  shows  that  overcasts  are 
usually  bluer  (i.e.,  have  higher  color  temperature) 
than  clear  daylight,6  although  this  enhanced  blue¬ 
ness  can  vary  significantly  from  minute  to  minute.  If 
our  claim  that  overcasts  are  bluer  than  clear  daylight 
seems  puzzling,  recall  that  Fig.  1  is  based  on  hemi¬ 
spheric  Eh  rather  than  on  narrow-FOV  Lx,  and  thus 
its  clear-sky  cases  are  dominated  by  direct  sunlight. 
During  the  day,  the  clear  sky’s  blue  is  indeed  much 
purer  than  that  of  the  bluest  overcast.  Third,  unlike 


Table  2.  Mean  Overcast  Chromaticities  and  Chromaticity  Gamuts 


Figure 

Dates 

Ex  or  Lxa 

h0  Interval  (°) 

Mean  u' 

Mean  v' 

i 

1 

9-5-02,  9-10-03 

Ex 

30.0-36.3 

0.20239 

0.47507 

0.0011029 

1 

11-8-03 

Ex 

30.0-36.3 

0.20291 

0.47238 

0.0076751 

1 

4-21-03 

Ex 

30.0-36.3 

0.19855 

0.46716 

0.015191 

3 

4-17-03 

u 

46.6-60.3 

0.20035 

0.46935 

0.0018261 

3 

5-22-03 

u 

46.6-60.3 

0.19707 

0.46780 

0.0013432 

3 

6-16-03 

u 

46.6-60.3 

0.19901 

0.46901 

0.0041030 

17 

10-12-02 

Ek 

5.0-12.3 

0.19893 

0.47027 

0.014782 

17 

11-15-03 

Ek 

5.0-12.3 

0.20014 

0.46677 

0.0070159 

19 

2-6-03 

u 

14.4-20.5 

0.19757 

0.46132 

0.0086214 

19 

2-17-03 

u 

14.4-20.5 

0.19824 

0.46523 

0.015537 

20 

1-1-03 

u 

12.3-20.2 

0.19690 

0.46519 

0.013589 

20 

1-2-03 

u 

12.3-20.2 

0.19709 

0.46276 

0.0066675 

20 

1-16-03 

u 

12.3-20.2 

0.19868 

0.46423 

0.0082847 

"Denotes  whether  chromaticities  in  a  given  figure  are  calculated  from  spectral  horizontal  irradiances  Ek  or  zenith  radiances  Lk. 
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Fig.  2.  Spectrally  integrated  horizontal  irradiances  E  versus  h0 
for  Fig.  l’s  two  overcasts.  Integrals  are  calculated  from  380-780  nm 
in  4  nm  steps. 


Fig.  4.  Spectrally  integrated  zenith  radiances  L  versus  h0  for  Fig. 
3’s  three  stratus  overcasts  at  Owings.  The  spectroradiometer’s 
measurement  FOV  is  1°,  and  integrals  are  calculated  from 
380-780  nm  in  4  nm  steps. 


the  often  smooth  chromaticity  curves  produced  by 
clear  twilights,17  daytime  overcast  curves  tend  to  be 
erratic,  with  colors  shifting  from  one  measurement  to 
the  next  in  response  to  clouds’  fluctuating  Ek  and 
optical  depth  t  (see  Fig.  2’s  time  series  of  integrated 
E).  Although  such  daylight  color  shifts  are  detectable 
in  principle,18  in  practice  most  go  unnoticed. 

Even  if  cloud  type  does  not  change,  the  color  gam¬ 
uts  of  overcasts  on  different  days  may  differ  percep¬ 
tibly  for  the  same  h0  interval.  Figure  3  shows 
chromaticities  from  zenith  Lx  for  three  different  stra¬ 
tus  overcasts  at  the  Owings  site  when  46.6°  <  h0 
<  60.3°.  Because  measured  overcast  chromaticities 
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Fig.  3.  Scatterplot  of  overcast  u' ,  v'  chromaticities  calculated 
from  zenith  spectral  radiances  Lx  at  Owings,  Maryland,  during 
stratus  overcasts  on  4-17-03,  5-22-03,  and  6-16-03;  h0  ranges  from 
46.6°-60.3°.  The  length  of  the  bar  labeled  “1  JND”  equals  the  local 
MacAdam  just-noticeable  difference  in  chromaticity. 


seldom  change  linearly,  beginning  in  Fig.  3  we  usu¬ 
ally  omit  line  segments  between  temporally  adjacent 
colors.  This  plotting  convention  makes  the  details  of 
our  chromaticity  gamuts  somewhat  clearer.  Like  all 
our  daytime  overcast  measurements  (i.e.,  those  taken 
when  h0  >  5°),  Fig.  3’s  chromaticities  are  to  the  left  of 
the  Planckian  locus.  Moreover,  these  three  overcasts 
have  more-or-less  distinct  chromaticity  clusters  that 
differ  from  one  another  on  the  CIE  diagram  by  ~1 
MacAdam  just-noticeable  difference  (or  JND;  1  JND 
is  the  horizontal  line  at  upper  right  in  Fig.  3). 19  In 
other  words,  if  typical  observers  simultaneously  com¬ 
pared  Fig.  3’s  clusters  of  overcast  colors,  most  could 
just  distinguish  among  them.  But  because  we  do  not 
actually  have  this  perceptual  luxury,  we  seldom  no¬ 
tice  such  small  color  shifts  between  different  over¬ 
casts. 

Figure  4,  which  shows  L  time  series  for  Fig.  3’s 
chromaticities,  helps  to  explain  differences  in  over¬ 
cast  color.  Although  all  three  overcasts  in  Figs.  3  and 
4  consisted  of  stratus,  their  appearance  ranged  from 
distinctly  variegated  (4-17-03)  to  almost  featureless 
(5-22-03).  As  clouds  moved  across  the  radiometer’s 
zenith  FOV,  the  resulting  changes  in  t  caused  the 
fluctuations  in  both  spectral  L,,  and  integrated  L  seen 
in  Figs.  3  and  4.  Because  the  5-22-03  overcast  had  the 
most  spatially  uniform  t,  it  has  the  smallest  L  dy¬ 
namic  range  (or  ratio  of  maximum  to  minimum  L ) 
and  smallest  g.  Conversely,  the  overcast  with  the 
largest  L  dynamic  range  (6-16-03;  see  Fig.  4’s  dashed 
curve)  has  the  largest^.  Thus  increased  variability  in 
overcast  t  tends  to  make  g  larger,  and  it  also  seems  to 
change  the  colorimetric  position  of  that  gamut.  In 
particular,  the  darkest  overcast  in  Fig.  4  (5-22-03)  is 
also  the  bluest  in  Fig.  3.  We  see  this  relationship 
directly  in  Fig.  5,  where  the  lower-E  (and  thus 
higher-T)  right  half  of  the  fisheye  image  is  distinctly 
bluer  than  the  higher-E  left  half.  Because  we  fixed 
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Fig.  5.  Composite  fisheye  image  of  a  bright  stratus  overcast  (left 
half,  E  =  24.3  W  rrT2)  and  a  much  darker  stratocumulus  overcast 
(right  half,  E  ~  3.6  W  nT2)  at  USNA  on  11-18-03  and  11-19-03, 
respectively.  Note  that  the  darker  overcast  is  distinctly  bluer.  Al¬ 
though  exposures  differ,  in  both  original  photographs  the  digital 
camera’s  white-balance  setting  was  the  same  and  h0  ~  18.5°. 


Fig.  6.  Correlated  color  temperature  (CCT)  in  Kelvins  versus  h0 
for  all  4278  of  our  horizontal  Ex  spectra  for  overcasts;  individual 
CCTs  are  plotted  as  dots.  The  right-hand  ordinate  gives  equivalent 
values  of  inverse  CCT,  which  is  measured  in  inverse  mega-Kelvins 
(106/CCT;  inverse  CCT  unit  is  MK1).  The  mean  CCT  curve  is 
calculated  using  bins  that  are  2°  wide  in  h0,  and  each  error  bar 
spans  2  standard  deviations  a  at  the  given  h0. 


the  digital  camera’s  white  balance  when  we  took  Fig. 
5’s  photographs,  their  color  shift  is  real  and  not  an 
electronic  artifact.  Optical  explanations  for  these 
overcast  color  phenomena  appear  in  Section  4. 

3.  Correlated  Color  Temperatures  of  Daytime 
Overcasts 

Correlated  color  temperature  (CCT)  has  long  been 
used  as  an  approximate,  convenient  alternative  to 
chromaticity  coordinates  near  the  Planckian  locus.20 
Because  CCT  reduces  two-dimensional  color  data 
(; u ',  a')  to  one  dimension,  CCT’s  convenience  comes  at 
the  cost  of  some  colorimetric  ambiguity.  For  example, 
although  a  6000-Kelvin  CCT  can  be  on  either  side  of 
the  Planckian  locus,21  all  our  measurements  of  day¬ 
time  overcast  CCTs  place  them  on  its  left  (or  green¬ 
ish)  side.  In  Fig.  6,  we  show  how  overcast  CCT  varies 
with  h0  for  nearly  half  our  9100  daytime  spectra. 
These  are  the  E x  spectra  measured  from  the  entire 
sky  hemisphere;  i.e.,  they  are  spectra  of  overcast  day¬ 
light  as  opposed  to  skylight.  Although  cloud  spectral 
optical  depth  tx  and  the  underlying  surface’s  spectral 
reflectance  rh  also  influence  overcast  CCT,4  h0  partly 
determines  the  daylight  spectrum  incident  on  an 
overcast’s  top.  In  turn,  this  daylight  spectrum  affects 
overcast  color  as  seen  from  the  ground. 

Because  the  incident  illumination  is  only  one  of 
several  factors  governing  overcast  color,  h0  and  CCT 
are  weakly  correlated  in  Fig.  6.  Yet  as  in  our  earlier 
work,6  mean  overcast  CCT  in  Fig.  6  clearly  has  a  local 
minimum  (i.e.,  is  reddest)  when  30°  <  h0  <  35°.  Al¬ 
though  Fig.  6’s  CCT  variance  might  eliminate  this 
local  minimum,  the  distinct  trend  in  mean  CCT  at 
nearby  h0  suggests  otherwise.  If  indeed  this  CCT 
minimum  is  real,  we  lack  a  ready  explanation  for  it: 


our  measurements  indicate  that  clear  daylight  at  the 
surface  (i.e.,  combined  skylight  and  direct  sunlight) 
usually  does  not  begin  to  redden  until  20°  <  h0 
<  25°.  Furthermore,  we  cannot  blame  the  overcast 
itself  for  this  reddening.  As  demonstrated  in  Section 
4,  diffuse  transmission  through  optically  thick  clouds 
makes  the  transmitted  light  slightly  bluer,  not  red¬ 
der. 

The  mean  and  median  overcast  CCTs  in  Fig.  6  are 
6358  K  and  6341  K,  respectively.  Although  these  val¬ 
ues  are  —375  K  larger  than  the  mean  overcast  CCTs 
that  we  reported  earlier,6  this  is  unlikely  to  be  per¬ 
ceptually  significant,  especially  given  that  the  overall 
standard  deviation  a  =  336  K  is  nearly  as  large  as 
the  CCT  difference.  In  fact,  our  new  mean  and  me¬ 
dian  CCT  are  closer  to  the  6500  K  overcast  white- 
balance  setting  suggested  by  some  digital  camera 
manufacturers.22  The  extreme  daytime  CCTs  are 
9316  K  (bluest)  and  5799  K  (least  blue),  and  we  plot 
the  corresponding  EK  spectra  in  Fig.  7.  The  high-CCT 
spectrum  occurred  when  h0  =  5.0°,  our  self-imposed 
lower  limit  on  daytime  h0.  That  the  maximum  CCT 
occurred  at  the  lowest  h0  is  not  surprising,  because 
even  though  direct  sunlight  atop  the  overcast  is  red¬ 
dened  then,  the  combined  effects  of  atmospheric  ex¬ 
tinction  and  the  cosine  law  mean  that  it  contributes 
relatively  little  to  horizontal  EK.17  Instead,  bluish  sky¬ 
light  dominates  illumination  both  above  and  below 
the  clouds. 

Figure  8’s  histogram  gives  the  distribution  of  all 
daytime  overcast  color  temperatures,  here  shown  as 
inverse  CCTs.6  The  unit  of  inverse  CCT  is  the  recip¬ 
rocal  mega-Kelvin  (106/CCT),  denoted  by  the  symbol 
MK  1  (originally  called  the  “mired”).23  Reciprocal 
mega-Kelvins  produce  a  uniform  scale  that  better 
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Fig.  7.  Normalized  horizontal  irradiances  Ek  for  the  two  daytime 
spectra  having  the  maximum  and  minimum  CCTs  in  our  data 
(9316  K  and  5799  K).  These  spectra  yield  our  colorimetric  ex¬ 
tremes  in  daytime  overcasts.  Each  original  spectrum  is  normalized 
by  its  sum;  to  convert  normalized  to  absolute  irradiances,  multiply 
the  11-22-03  Ek  by  0.814  W  rrT2  nrrT1  and  the  11-8-03  Ek  by 
16.3  W  rrT2  nm  \ 


Fig.  9.  Mean  T(\)  transmission  spectra  for  Ek  spectra  measured 
at  USNA  during  stratus  overcasts  on  4-4,  4-8,  and  4-21-03;  h0 
ranges  from  40.1°-57.7°.  N  =  25  and  N  =  17  irradiance  spectra 
were  used  to  calculate  the  mean  T(X,  low)  and  T(\,  high)  spectra, 
respectively.  Each  error  bar  spans  2<r  at  the  given  wavelength. 
Best-fit  equations  are  for  the  interval  400  £  X  <  680  nm,  with 
proportionality  constants  kk  =  4.077  and  k2  =  1.946. 


describes  human  color  sensitivity  than  does  CCT  it¬ 
self.  On  this  inverse  CCT  scale,  our  overcast  spectra 
have  a  mean  =  157.7  MKT1,  ct  =  7.79  MK_1,  and  a 
modal  interval  of  154-156  MKT1.  Clearly  Fig.  8  is 
skewed  left  toward  smaller  inverse  CCT,  and  this  is 
caused  by  a  relatively  few  high-CCT  overcasts  seen  at 
low  h0  (see  Fig.  6’s  left  side). 

4.  Blue  Clouds  and  Overcast  Transmission  Spectra 

A  fundamental  problem  in  analyzing  overcast  colors 
is  that  the  continual  changes  in  cloud  E  and  L  (see 
Figs.  2,  4)  further  complicate  an  already  complex  op- 


Fig.  8.  Histogram  of  inverse  CCT  for  our  4278  overcast  Ex  spec¬ 
tra.  The  modal  interval  is  154-156  MKT1,  and  each  bin  is  2  MET1 
wide. 


tical  system.  We  take  two  different  tacks  in  address¬ 
ing  this  problem:  (1)  averaging  measured  spectra 
over  a  small  range  of  h0  and  (2)  modeling  the  effects 
of  individual  parameters  on  overcast  colors. 

We  start  by  defining  an  overcast’s  spectral  trans¬ 
missivity  T(\),  where 


E(X,  OVC) 
T(X)  =  E(X,  CLR) 


(1) 


is  the  ratio  of  an  overcast’s  downwelling  spectral  ir¬ 
radiances  E(X,  OVC)  to  those  of  some  average 
E{X,  CLR)  for  clear  daylight  at  the  same  h0.  This  last 
qualification  is  needed  because  both  E(X,  OVC)  and 
E(X,  CLR)  are  functions  of  h0.  Here  we  average  Eh 
spectra  from  the  clear  afternoons  of  9-5-02  and 
9-10-03  at  USNA  (see  Fig.  1).  In  radiative-transfer 
models,  an  overcast’s  T  usually  is  usually  defined  as 
the  ratio  of  downwelling  irradiance  at  cloud  bottom  to 
that  at  cloud  top  (i.e.,  incident  daylight).  However, 
our  ground-based  measurement  of  T  is  nearly  equiv¬ 
alent  because  many  overcast  cloud  tops  and  bottoms 
are  fairly  close  to  the  ground.24 

Figure  9  shows  two  stratus  T(X)  spectra,  each  av¬ 
eraged  over  N  overcast  irradiance  spectra  with  either 
relatively  low  T  (N  =  25)  or  high  T  (N  =  17).  The 
original  overcast  Ex  spectra  were  acquired  over  a  re¬ 
stricted  h0  range  on  several  different  days.  To  show 
more  spectral  detail,  we  draw  the  two  T{X)  spectra 
with  different  ordinate  scales,  and  these  scales  pre¬ 
serve  the  relative  slopes  of  the  two  curves.  Figure  9 
provides  unambiguous  spectral  evidence  that  both 
low-71  (or  darker)  and  high-71  (or  brighter)  overcasts 
shift  incident  daylight  toward  the  blue,  as  indicated 
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Fig.  10.  Mean  T(k)  spectra  for  Ex  spectra  measured  at  USNA 
during  a  stratocumulus  overcast  on  3-19-03;  h0  ranges  from 
41.9°-50.5°.  N  =  15  and  N  =  19  irradiance  spectra  were  used  to 
calculate  the  mean  low-  and  high-transmissivity  T(k)  spectra,  re¬ 
spectively.  Each  error  bar  spans  2<r  at  the  given  wavelength,  and 
for  clarity  we  show  bars  only  for  the  T(k,  high)  curve.  Best-fit 
equations  are  for  the  interval  400  <  X  £  680  nm,  with  proportion¬ 
ality  constants  k3  =  4.136  and  kA  =  3.943. 


by  the  best-fit  power-law  relationships  T(k)  °c  \  0  611 
for  low  T  and  T(k)  oc  X  0  27(1  for  high  T.  A  slightly  more 
pronounced  X  '  dependence  exists  for  stratocumulus 
overcasts  (see  Fig.  10).  For  example,  at  comparable  T 
values  in  Figs.  9  and  10,  c  is  —44%  larger  for  the 
stratocumulus  overcasts.  Thus  at  the  same  T  or  t,  a 
stratocumulus  overcast  will  look  slightly  bluer  than  a 
stratus  overcast.  However,  this  distinction  is  often 
obscured  for  two  reasons:  (1)  T  is  seldom  equal  in 
different  overcasts  (see  Fig.  2),  and  (2)  temporal  color 
constancy25  makes  us  less  sensitive  to  any  such  daily 
changes  in  overcast  color. 

At  wavelengths  >550  nm,  the  measured  T(X)  spec¬ 
tra  in  Figs.  9  and  10  have  local  maxima  and  minima 
similar  to  those  found  in  transmission  spectra  asso¬ 
ciated  with  continuum  absorption  by  water  vapor  and 
oxygen.26  Thus  absorption  by  atmospheric  gasses 
within  clouds  adds  some  spectral  detail  to  overcast 
T(X)  spectra,  but  it  does  not  explain  their  X  '  shape. 

The  fact  that  overcast  T{\)  varies  approximately  as 
k~c  may  be  surprising,  since  this  is  also  true  of  the 
Rayleigh  blue  sky,  for  which  Lx  A.  4  (i.e.,  c  =  4).  Yet 
unlike  the  clear  sky,  in  an  overcast  the  amount  of 
bluing  increases  with  t  (i.e.,  with  decreasing  T).  Thus 
thicker  clouds  are  bluer.  This  strongly  suggests  that 
spectrally  selective  absorption  rather  than  scattering 
ultimately  causes  the  subtle  bluing  seen  beneath 
overcasts.27  In  fact,  we  can  use  T  to  estimate  overcast 
t  by  modifying  Bohren’s  two-stream  radiative- 
transfer  model28  to  include  surface  reflectance  r. 


wavelength  X  (nm) 

Fig.  11.  MODTRAN4  simulated  T(k)  spectra  for  the  model’s  de¬ 
fault  stratus  overcast  when  its  clouds  are  100  m  or  670  m  thick.  In 
both  cases,  h0  =  45°,  and  cloud  base  is  330  m  above  the  surface.  For 
the  thicker  overcast,  we  calculate  T(k)  using  combined  direct-beam 
and  diffuse  irradiances  at  the  surface,  but  we  switch  to  irradiances 
calculated  at  cloud  top  and  bottom  for  the  thinner  overcast  (see  Eq. 
(1)).  Best-fit  equations  are  for  the  interval  380  £  \  <  780  nm,  with 
proportionality  constants  k$  =  1.069  and  ke  =  2.544. 

Then 

2C r1  - 1) 


where  g  is  the  single-scattering  asymmetry  parame¬ 
ter  for  cloud  droplets,  and  all  variables  are  implicit 
functions  of  X.  Not  surprisingly,  Eq.  (2)  is  most  real¬ 
istic  for  thicker  overcasts  where  the  assumptions  of 
the  two-stream  model  are  best  satisfied. 

5.  Modeling  Overcast  Transmission  Spectra  and 
Chromaticities 

To  analyze  the  fundamental  optical  properties  that 
determine  overcast  T{\),  we  use  the  detailed  and 
extensively  tested  radiative-transfer  model  MOD- 
TRAN4.  MODTRAN  combines  features  of  horizon¬ 
tally  homogeneous  spherical-shell  and  plane-parallel 
atmospheres,  and  it  calculates  solar  direct-beam  and 
diffusely  scattered  components  at  each  of  the  model’s 
many  vertical  levels.  Provided  that  h0  >  0°,  MOD¬ 
TRAN  yields  quite  realistic  Lx  and  Ex  spectra  for 
skylight  and  daylight.29  Except  as  noted  below,  each 
of  our  MODTRAN  simulations  includes  the  following 
parameters:  a  default  midlatitude  winter  atmo¬ 
spheric  profile  of  temperature,  pressure,  humidity, 
and  gas  mixing  ratios;  tropospheric  aerosols  typical  of 
rural  sites;  background  stratospheric  dust  and  other 
aerosols;  multiple  scattering;  a  Lambertian  surface 
with  r  =  0.2  at  all  X;  and  Mie  aerosol  phase  functions. 

In  Fig.  11,  we  calculate  MODTRAN’ s  T(k)  using 
downwelling  surface  Eh  for  a  stratus  layer  whose 
thickness  A z  =  670  m  and  clear-sky  Ex  at  the  same 
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Fig.  12.  Comparison  of  simulated  MODTRAN  T(\)  for  the  direct- 
beam  and  diffuse  components  in  Fig.  ll’s  100  m  thick  stratus 
overcast.  T(\,  diffuse)  >  1  because  the  diffuse  irradiances  at  cloud 
bottom  (the  numerator  in  Eq.  (1))  include  energy  scattered  from 
the  direct  solar  beam,  whereas  diffuse  irradiances  at  cloud  top 
exclude  direct  sunlight. 


Fig.  13.  MODTRAN  visible-wavelength  spectra  of  absorption 
cross  section  Cabs  and  extinction  cross  section  Cext  for  stratus  drop¬ 
lets  of  radius  a.  Cext  is  normalized  by  its  value  at  550  nm  and  is 
proportional  to  2ira2.  Note  that  Cabs  is  scaled  logarithmically. 


h0.  We  also  include  in  Fig.  11a  stratus  overcast  with 
A z  =  100  m,  but  here  we  calculate  T{X)  using  simu¬ 
lated  Ex  at  cloud  top  and  bottom.  This  change  lets  us 
examine  in  Fig.  12  the  diffuse  and  direct-beam  T{X) 
restricted  to  just  the  cloud  itself.  For  the  100  m  over¬ 
cast  in  Fig.  11,  mean  T  is  naturally  rather  large,  and 
T(X)  is  nearly  spectrally  neutral  (c  =  0.065).  Yet  un¬ 
like  many  other  models,  MODTRAN’s  simulated  T(X) 
spectra  closely  resemble  our  measured  T(X)  spectra, 
including  such  features  as  continuum  absorption  by 
water  vapor  and  oxygen.  For  example,  compare  the 
similar  transmissivity  curves  T{ X,  670  m)  in  Fig.  11 
and  T(X,  high)  in  Fig.  9. 

Figure  12  shows  the  very  different  MODTRAN  T(X) 
spectra  of  daylight’s  attenuated  direct-beam  and  am¬ 
plified  diffuse  scattering  components.  Even  a  thin 
stratus  overcast  with  A z  =  100  m  nearly  extinguishes 
the  direct  solar  beam  and  makes  it  slightly  bluer,  as 
seen  in  the  T(X,  direct)  spectrum.  Yet  Fig.  12’s 
T{X,  diffuse)  spectrum  is  both  >1  and  redder  at  the 
cloud’s  bottom.  This  occurs  because,  by  definition,  the 
model’s  cloud-top  diffuse  Ex  consists  only  of  energy 
from  outside  the  Sun’s  solid  angle.  Thus  it  starts  with 
a  diffuse  EK  component  predominated  by  blue  sky¬ 
light,  but  after  traversing  the  cloud,  this  diffuse  Ex 
has  ( 1)  gained  energy  scattered  from  the  direct  beam, 
which  (2)  was  much  less  blue  than  skylight  to  begin 
with,  and  from  which  (3)  there  is  slightly  preferential 
scattering  of  red  light  that  adds  to  the  diffuse  down¬ 
ward  Ek. 

In  Fig.  13,  we  examine  how  absorption  and  extinc¬ 
tion  by  stratus  cloud  droplets  can  account  for  MOD¬ 
TRAN’s  overcast  colors.  First,  note  that  the  size- 
dependent  extinction  for  these  relatively  large 
(radius  a  >  5  p,m)  droplets  is  simply  proportional  to 
2ua2.  Second,  MODTRAN’s  absorption  coefficient  for 
liquid  water  (Cabs)  is  significant  only  at  X.  >  590  nm, 


which  apparently  leads  to  droplets  that  have  slightly 
nonneutral  extinction  cross  sections  Cext  at  visible  X 
(Cext  increases  ~2.7%  from  380-780  nm). 

So  if  extinction  increases  weakly  with  X,  then  the 
(possibly)  transmitted  direct  beam  becomes  slightly 
bluer,  and  the  scattered  light  becomes  slightly  less 
blue,  as  we  saw  earlier.  How  do  these  droplet-level 
details  help  explain  thick  clouds’  bluish  color?  Imag¬ 
ine  an  optically  thick  overcast  with  cloud  normal  t 
>  100.  We  routinely  experience  such  overcasts,  and 
as  seen  in  Fig.  14,  surface  E  beneath  them  decreases 
monotonically  (but  not  linearly)  with  increasing  A z  or 


Fig.  14.  MODTRAN  simulation  of  horizontal  irradiances  E  at  the 
surface  as  a  function  of  cloud  thickness  A z.  For  normal  optical 
depth  t,  the  E( t)  curve  is  not  exactly  congruent  with  the  E(Az) 
curve  because  we  calculate  the  former  independent  of  MODTRAN 
using  Eqs.  (1)  and  (2). 
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Fig.  15.  MODTRAN  simulations  of  stratus  overcast  u' ,  v'  chro- 
maticities  calculated  from  Ex  as  functions  of  A z  and  h0.  The  dashed 
h0  curve  sets  Az  =  1  km,  and  the  solid  Az  curve  sets  h0  =  45°.  The 
underlying  Lambertian  surface  has  a  constant  spectral  reflectance 
rx  =  0.2.  During  the  daytime,  simulated  overcast  colors  grow 
steadily  bluer  as  Az  increases  and  h0  decreases. 

t.  In  Fig.  14,  we  use  Eqs.  (1)  and  (2)  to  calculate  t  from 
MODTRAN  overcast  E. 

Clearly  in  a  neutrally  absorbing,  optically  thick 
cloud,  multiple  scattering  by  large  drops  (i.e.,  a  k) 
would  only  yield  TCk)  =  constant.  As  in  the  popular 
misconception,  light  exiting  such  a  cloud’s  bottom 
would  have  the  same  color  as  the  daylight  incident  on 
its  top.  However,  in  real  clouds  with  nonneutral  Cabs, 
multiple  scattering  for  large  Az  or  t  amplifies  spec¬ 
trally  selective  absorption  enough  to  produce  the 
clouds’  bluish  colors.27  Because  this  bluing  increases 
with  cloud  t,  absorption  within  water  droplets  and 
not  droplet  scattering  in  itself  likely  causes  overcast 
colors.  Multiple  scattering  also  increases  absorption 
by  the  increased  amount  of  water  vapor  in  and  near 
the  clouds,  although  this  does  not  explain  the  overall 
\-c  shape  of  T(X). 

Now  we  use  MODTRAN  to  help  understand  how 
overcast  colors  depend  on  h0  and  cloud  Az.  In  Fig.  15, 
MODTRAN’s  simulated  stratus  overcast  colors  be¬ 
come  progressively  bluer  as  Az  increases  from  0  to 
10  km,  with  CCT  increasing  by  —1660  K  for  h0 
=  45°.  For  fixed  Az  =  1  km,  an  even  larger  color  shift 
occurs  as  h0  decreases  from  45°  to  5°  in  Fig.  15.  Unlike 
Fig.  6’s  observed  local  minimum  in  CCT  near  h0 
—  30°,  MODTRAN  predicts  only  a  steady  increase  in 
CCT,  with  much  of  that  increase  occurring  for  5° 
<  h0  <  10°.  In  real  overcasts,  however,  CCT  can 
actually  decrease  with  decreasing  h0  if  Az  also  de¬ 
creases  then  (e.g.,  see  the  h0  extremes  in  Fig.  l’s 
Granada  11-8-03  chromaticity  curve). 

Might  thick  overcasts  be  bluish  because  their  tops 
are  illuminated  by  clear  skylight  whose  blueness  in¬ 
creases  with  altitude  z?  We  know  that  skylight  purity 
increases  with  z  because  multiple  scattering  in  the 
clear  atmosphere  decreases  with  z.30  Yet  both  a  sim- 


Fig.  16.  MODTRAN  simulations  of  stratus  overcast  u',v’  chro- 
maticities  calculated  from  Ek  as  functions  of  Az  for  the  cases  of 
fixed-altitude  tops  z(top)  and  fixed-altitude  bottoms  z(base).  Al¬ 
though  the  Case  2  chromaticities  (fixed,  high-altitude  z(top))  are 
slightly  bluer,  in  both  cases  bluing  of  daylight  beneath  the  overcast 
depends  much  more  strongly  on  cloud  Az. 

pie  thought  experiment  and  MODTRAN  simulations 
show  that  cloud-top  height  z(top)  contributes  little  to 
overcast  bluing.  Certainly  downwelling  clear-sky  Ex 
grow  bluer  as  observer  z  increases,  but  only  within 
limits.  For  example,  we  have  measured  daytime  Ek 
spectra  at  altitudes  of  z  —  3.8  km  that  are  distinctly 
bluer  than  those  near  sea  level.  However,  this  trend 
cannot  continue  indefinitely,  since  skylight  radiances 
decrease  with  z  even  as  their  colorimetric  purities 
increase.  Ultimately,  the  purest  blue  skylight  occurs 
at  z  such  that  only  one  molecule  scatters  sunlight  to 
us;  but  then  the  sky  will  be  black,  not  blue,  and  the 
color  of  daylight  will  be  that  of  extraterrestrial  sun¬ 
light. 

In  fact,  blue  skylight’s  increasing  purity  with  z  has 
only  a  minor  effect  on  overcast  colors  at  the  surface. 
Depending  on  z(top),  daylight  incident  on  higher 
clouds  may  well  be  bluer  than  that  on  lower  clouds, 
but  the  clear-sky  molecular  atmosphere  still  exists 
between  the  surface  and  z(top).  Although  the  illumi¬ 
nation  geometry  certainly  changes  for  gas  molecules 
within  clouds,  their  spectral  scattering  properties  are 
unchanged,  and  so  will  be  their  individual  contribu¬ 
tions  to  surface  Ek.  MODTRAN  simulations  show 
that  cloud  thickness  Az  plays  a  far  more  important 
role  than  z(top).  In  Fig.  16,  we  plot  as  functions  of  Az 
the  daylight  chromaticities  beneath  two  MODTRAN 
stratus  overcasts  for  h0  =  45°.  These  overcasts  have 
either  a  fixed  bottom  z(base)  =  0.33  km  (Case  1)  or  a 
fixed  z(top)  =  5.33  km  (Case  2);  all  other  MODTRAN 
parameters  are  identical.  As  cloud  Az  increases  from 
right  to  left  in  Fig.  16,  z(top)  either  increases  (Case  1) 
or  is  constant  (Case  2).  Despite  distinct  differences  in 
daylight  spectra  atop  the  Case  1  and  Case  2  over¬ 
casts,  they  yield  nearly  identical  chromaticity  curves. 

If  daylight  color  beneath  these  two  overcasts  were 
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Fig.  17.  Temporal  trends  in  overcast  u',v'  chromaticities  calcu¬ 
lated  from  Ex  near  Marion  Center,  Pennsylvania,  on  10-12-02  and 
at  Granada  on  11-15-03;  h0  ranges  from  5.0°-12.3°. 
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Fig.  18.  MODTRAN  simulations  of  stratus  overcast  u',v'  chro¬ 
maticities  at  the  zenith  as  functions  of  cloud  A z  and  the  underlying 
surface’s  spectral  reflectance  rx.  The  thick  arrow  indicates  MOD- 
TRAN’s  prediction  of  the  chromaticity  shift  caused  by  changing  the 
surface  from  grass  to  snow  for  Fig.  19’s  observed  changes  in  cloud 
A z.  For  historical  purposes,  we  also  include  overcast  chromaticities 
simulated  by  Middleton.4 


largely  determined  by  z(top),  then  chromaticity  dis¬ 
tances  Au'v’  at  the  same  A z  should  be  nearly  as  large 
as  or  larger  than  those  caused  by  increasing  cloud  Az 
(pairs  of  +’s  in  Fig.  16  mark  equal-Az  chromaticities). 
Yet  Fig.  16  clearly  shows  that  this  is  not  so:  although 
Case  2  chromaticities  are  slightly  bluer  than  those  for 
Case  1  at  the  same  Az  (i.e.,  Case  2  chromaticities  are 
shifted  slightly  leftward),  the  maximum  Au'v'  at  Az 
=  2.5  km  is  far  less  than  the  Au'v'  caused  by  a  2.5  km 
thick  stratus  overcast  (i.e.,  the  chromaticity  distance 
between  Az  =  2.5  km  and  Az  =  10  m).  In  particular, 
note  that  the  largest  Case  1-Case  2  color  difference  at 
Az  =  2.5  is  still  less  than  1 JND.  Similar  results  occur 
at  smaller  h0.  Thus  multiple  scattering  and  absorp¬ 
tion  within  clouds,  not  the  altitudes  of  cloud  tops, 
chiefly  determine  daytime  overcast  colors  at  the  sur¬ 
face. 

In  one  sense,  Fig.  17’s  measured  chromaticities 
support  MODTRAN’s  claim  in  Fig.  15  about  the 
/in-dependence  of  overcast  colors:  for  both  overcasts, 
the  largest  CCT  increases  occur  at  the  smallest  h0.  In 
Fig.  17,  CCT  increases  steadily  only  when  h0  <  8.2° 
(Granada  11-15-03)  or  h0  <  6.6°  (Marion  Center  10- 
12-02).  However,  the  measured  effects  of  cloud  Az 
seem  to  differ  from  those  predicted  by  MODTRAN. 
Although  the  Granada  overcast  is  consistently  much 
darker  than  that  at  Marion  Center  (by  an  E  factor 
— 2. 1),  its  mean  u'  is  in  fact  larger,  not  smaller  as 
MODTRAN  predicts.  Furthermore,  MODTRAN  pre¬ 
dicts  that  the  Granada  chromaticity  curve  will  shift 
rightward  <scl  JND  as  a  result  of  Granada  having  a 
stratocumulus,  rather  than  a  stratus,  overcast.  In¬ 
stead,  Fig.  17’s  Granada  curve  >1  JND  to  the  right  of 
the  Marion  Center  curve.  Our  point  here  is  not  that 
MODTRAN  is  unrealistic  but  that  ultimately  we 
must  know  more  about  overcast  illumination  and  in¬ 
ternal  structure  (e.g.,  clouds  at  other  z,  vertical  inho¬ 
mogeneities  in  droplet  and  haze  number  densities, 


detailed  drop  size  distributions)  before  we  can  test 
the  model  rigorously. 

6.  Overcast  Color  and  Surface  Color 

Some  50  years  ago,  Middleton  made  one  of  the  first 
systematic  attempts  to  predict  overcast  sky  colors.  He 
did  so  by  incorporating  spectral  ground  reflectance  rx 
and  scattering  and  absorption  by  cloud  droplets  into 
a  model  of  overcast  spectra  in  a  plane-parallel  atmo¬ 
sphere.  Other  model  parameters  included  cloud 
thickness,  droplet  number  density,  and  monodisperse 
droplet  radius  and  scattering  coefficient;  h0  was  fixed 
at  23°.  When  Middleton  ultimately  calculated  chro¬ 
maticities  from  Ek  beneath  the  overcast,  he  decided  to 
ignore  absorption  by  both  droplets  and  their  inclu¬ 
sions,  asserting  that  the  overcast  sky’s  chromaticity 
and  luminance  are  to  “a  considerable  extent  a  func¬ 
tion  of  the  reflectance  of  the  ground,  but  depend  only 
to  a  minor  degree  on  the  presence  of  dissolved  ab¬ 
sorbing  matter  in  the  cloud  water.”4  As  beneficiaries 
of  several  additional  decades  of  research  on  overcasts, 
what  new  can  we  say  about  Middleton’s  claims? 

In  Fig.  18,  the  short  line  segment  at  upper  right 
connects  Middleton’s  predicted  chromaticities  for 
overcasts  with  cloud  Az  =  1  km  or  2  km  above  grass. 
A  longer  line  at  Fig.  18’s  bottom  connects  his  chro¬ 
maticities  for  the  same  Az  pair  above  snow.  Middle- 
ton  thus  predicts  that  the  same  overcast  will  be 
markedly  bluer  above  snow  than  above  grass.  Al¬ 
though  MODTRAN  makes  the  same  qualitative  pre¬ 
diction  about  snow’s  effect,  quantitatively  it  greatly 
expands  Middleton’s  overcast  g  for  overcast/grass 
and  shifts  those  chromaticities  toward  higher  CCT. 
Given  that  both  models  use  similar  rx  spectra  for 
grass  and  for  snow,  their  differences  probably  are  not 
caused  by  the  underlying  surfaces. 

Because  snow  reflects  more  bluish  cloud  light  to  the 
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Fig.  19.  Scatterplot  of  overcast  u' ,  v'  chromaticities  calculated 
from  zenith  Lk  at  Owings  during  stratus  overcasts  on  2-6-03  and 
2-17-03;  h0  ranges  from  14.4°-20.5°.  Only  on  2-17-03  was  the 
ground  snow-covered.  Compare  these  measured  chromaticities 
with  MODTRAN  simulations  in  Fig.  18. 


cloud’s  bottom  than  does  grass  (i.e.,  the  cloud  bottom 
is  brighter  above  snow),  there  is  more  scattering 
within  the  cloud.  More  germane  for  us  is  that  the 
additional  scattering  results  in  additional  spectrally 
selective  absorption  by  the  cloud.  Some  of  the  remain¬ 
ing  blue-enhanced  light  is  then  scattered  downward 
to  our  eyes  and  radiometers,  making  the  overcast 
both  bluer  and  brighter  than  it  would  be  over  a  low-r 
surface.  In  effect,  snow  increases  overcast  blueness 
by  increasing  cloud  Az  and  t. 

Middleton’s  chromaticities  for  an  overcast  above 
snow  clearly  are  incorrect  because  they  are  on  the 
opposite  side  of  the  Planckian  locus  from  our  mea¬ 
sured  chromaticities.  This  problem  may  be  caused  by 
Middleton’s  decision  to  ignore  spectral  absorption  by 
cloud  water.  Yet  MODTRAN  avoids  this  problem 
above  snow  cover  only  to  encounter  others.  In  Fig.  19, 
the  2-17-03  overcast  immediately  followed  a  major 
winter  storm  that  left  —0.3  m  of  snow  on  the  ground 
at  Owings;  on  2-6-03  the  ground  was  snow-free.  On 
both  days,  nearly  featureless  stratus  clouds  com¬ 
prised  the  overcast.  The  most  obvious  colorimetric 
differences  between  Fig.  19’s  two  overcasts  are  that 
the  post-snow  g  increases  in  size  and  shifts  slightly 
leftward.  This  observed  behavior  is  quite  different 
from  that  of  MODTRAN,31  which  predicts  large  shifts 
away  from  the  Planckian  locus  and  toward  higher 
CCTs  (see  Fig.  18’s  thick  arrow). 

In  fact,  surface  rA  may  matter  less  than  either 
Middleton’s  model  or  MODTRAN  suggest.  Using  the 
same  u',v'  scale  as  Figs.  18  and  19,  Fig.  20  shows 
overcast  zenith  chromaticities  measured  at  Owings 
on  three  snow-free  days  in  January  2003.  Although 
the  surface  rk  changed  little  during  this  two-week 
period,  the  spread  in  Fig.  20’s  overcast  chromaticities 
is  comparable  to  that  in  Fig.  19.  Furthermore,  the 
1-1-03  overcast  has  many  chromaticities  that  are  dis¬ 


Fig.  20.  Scatterplot  of  overcast  u' ,  v'  chromaticities  calculated 
from  zenith  Lk  at  Owings  during  stratus  overcasts  on  1-1-03,  1-2- 
03,  and  1-16-03.  The  ground  was  snow-free  on  all  days;  h0  ranges 
from  12.3°-20.2°. 


placed  even  farther  from  the  Planckian  locus  than  are 
Fig.  19’s  most  distant  post-snow  colors.  This  observa¬ 
tion  naturally  raises  the  question  of  how  much  of  Fig. 
19’s  chromaticity  shift  is  due  to  snow.  To  us,  Figs.  19 
and  20  suggest  that  daily  (and  even  minute-by- 
minute)  changes  in  cloud  tx  may  well  do  more  to  affect 
overcast  colors  than  all  but  the  most  extreme  shifts  in 
surface  rk. 

7.  Conclusions 

Our  measurements  of  overcast  T(\)  spectra  (Figs.  9 
and  10)  should  help  dispel  any  lingering  notion  that 
overcasts  are  spectrally  neutral  transmitters  of 
visible-wavelength  daylight.  While  overcast  colors  lit¬ 
erally  pale  in  comparison  to  their  spectacular  cer¬ 
ulean  counterparts,  they  have  much  to  teach  us  about 
radiative  transfer  in  clouds.  Although  the  total  depth 
of  liquid  water  in  even  the  thickest  overcast  is  quite 
small,  multiple  scattering  greatly  enhances  the  small 
amount  of  spectrally  selective  absorption  that  occurs 
each  time  light  interacts  with  a  cloud  droplet.28  After 
many  scatterings,  the  amount  of  bluing  caused  by 
thick  clouds  is  quite  perceptible,  as  even  a  casual 
inspection  of  their  bottoms  reveals  (Fig.  5).  Contin¬ 
uum  absorption  by  water  vapor  and  oxygen  adds  im¬ 
portant  spectral  details  at  longer  visible  wavelengths 
in  T(k)  spectra  (Figs.  9  and  10), 26  but  does  not  explain 
their  overall  shape  or  the  resulting  overcast  colors. 

However  drab  the  colors  of  overcast  skies  may 
seem  visually,  explaining  their  chromaticities  and 
(ir)radiance  spectra  continues  to  be  a  challenge.  First, 
the  observed  g  of  daylight  from  overcasts  usually  is 
larger  than  that  of  clear-sky  daylight  (Fig.  1),  and 
MODTRAN  simulations  only  begin  to  mimic  this 
range  of  colors  (Figs.  15  and  18).  Second,  sorting  out 
the  relative  spectral  contributions  of  overcasts’  T(\) 
or  tx,  the  Ex  spectra  incident  on  the  clouds’  tops,  and 
the  rx  spectra  of  the  underlying  surface  remains  dif- 
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ficult  many  years  after  Middleton’s  pioneering  work. 
Our  research  suggests  that  rk  probably  is  less  impor¬ 
tant  than  he  concluded,  and  that  other  purely  atmo¬ 
spheric  factors  such  as  the  presence  of  fog  (see  Fig. 
20’s  1-1-03  overcast)  are  more  important.  Finally, 
overcasts’  apparent  visual  simplicity  hides  great  col¬ 
orimetric  and  spectral  complexity.  To  understand 
this  complexity  better,  in  future  research  we  will  con¬ 
sider  vertical  and  horizontal  inhomogeneities  within 
optically  thick  clouds.  Thus  while  overcast  colors  may 
look  pedestrian,  they  offer  a  wealth  of  new  insights  in 
atmospheric  optics. 
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